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On . Abstract 

a^ ■ 

T-H ■ The influence of an external electric field on the spin dynamics of an electrically neutral Fermi 
liquid is considered, the mechanism of such an influence being the relativistic spin-orbital interac- 

^ I tion. As a result, Leggett's equations for the spin dynamics of weakly polarized Fermi liquids are 
generalized to the case of non-zero external electric field. In addition, we obtained the transverse 

00 ■ spin dynamics equation for strongly spin-polarized liquids in an electric field at zero temperature. 

', In both situations covariant derivatives depending on the electric field are shown to be substituted 

I for spatial gradients in line with the SU(2) gauge invariance of the microscopic Hamiltonian. The 

J> ■ new equations are applied to the study of spin flow along a channel, where an electric field is found 

, to bring about an additional phase shift of the order of magnitude of the phase shift in superfluid 

[ ^He-S but growing with time. 
00 ■ 

^ . PACS numbers: 67.55.-s, 67.65. -|-z 
0\ ' 



g 1 Introduction 
I 

^ The advances in experimental technique in recent years made feasible the observation 
Q , of effects due to electric field on the motion of electrically neutral superfluid ^He [|T| . 
> ! Experiments on superfluid ^He-5 in an electric field, using the spin current Josephson 
^ effect, are under way at the moment p. The influence of an external electric field 
on superfluid ^He has also been treated theoretically in scientific literature [0, 0] . 
However, a similar investigation for a normal Fermi liquid is lacking. It is this question 
which is attempted to be answered by the present paper. 

There have been predicted two competitive mechanisms of the influence of an elec- 
tric field on the spin dynamics of the superfluid ^He viz., (i) slight deformation of 
the electronic shells of ^He atoms by the gradient of the order parameter and (ii) spin- 
orbital coupling [fi X p)E /mc of the electric field with the moving magnetic moments 
of the ^He nuclei, where /Lt is the magnetic moment of ^He atom. In the present paper 
this latter relativistic effect is incorporated into spin dynamics equations of a normal 
Fermi liquid. The result complies with the gauge invariance arguments developed in 
Ref. 
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01 lengtn connecting two reservoirs witn nomogeneousiy precessmg magnetizations, 
the supercurrent, induced by an electric field, leads to an additional phase shift of the 
order of 

5a ~ - 10-4 rad (1) 

he 

for E" « 3 • 10'^ V/cm and L ^ 1 cm. According to this is an experimentally 
detectable quantity. 

In a normal Fermi liquid a phase shift of the same order is demonstrated to arise 
and grow with time in a similar one-dimensional spin fiow geometry. 

The set-up of the paper is as follows. The second section is dedicated to introducing 
the SU(2) gauge invariant theory of the interaction of electromagnetic field with a 
Fermi liquid. In the third section the weakly polarized Fermi liquid is considered on 
the basis of the modifications to the Landau theory resulting from the appearance of 
an electric field. Although this section does not overtly leans on the gauge invariance, 
the outcome agrees with what one expects from Sec. II. The next section is given 
over to the study of one-dimensional spin fiow as an example of application of the new 
equations. 

The invariance properties of the microscopic Hamiltonian are employed in treating 
strongly spin-polarized liquids at zero temperature in Sec. V. The last section sums 
up the main results. 



2 Gauge Invariance 

In the SU(2)-gauge- invariant wording of Ref. [jSj] electric and magnetic fields are intro 
duced uniformly as components of a gauge field A^, where /i = 0, 1, 2, 3: 

9 



gAo = gH = uji^, gA^ 



c 



(2) 



Here g is the gyromagnetic ratio for ^He nuclei, u^l is the Larmor frequency. 

The Hamiltonian of particles with the mass m and the magnetic moment /i = hg/2 
in an external electromagnetic field in the second quantization representation is 

fi(H - —p X E 



Pt 



2m 



mc 



(3) 



where fi = fi& is the magnetic moment operator. The Hamiltonian can be rewritten 
equivalently using the gauge field A^: 



2m 



^/^(r) d r. 



(4) 



The last expression is invariant against rotations of the spin space if only they are 
accompanied by the corresponding change of the gauge field: 



tjj n^Jj, gA^ g^A,, + a;„. 



(5) 
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u ana /-Cyu) is a oJJ rotation matrix corresponamg to /-c. ine rrequency tensor or tnis 
rotation u;^ is defined as = IZ^d^lZ^ and the gradients are designated as d^^ 

where Oq = dt, di = Vi. 

This symmetry allows one to refer to as the gauge field. Quite generally, gauge 
fields appear when the invariance of a Hamiltonian against some local gauge symme- 
try is stipulated. Nevertheless, it might be worth mentioning that this symmetry is 
somewhat formal since it is valid only for the Hamiltonian of ^He atoms and not for 
the electromagnetic field. 

The SU(2) gauge invariance of the Hamiltonian calls for all observables to be gauge 
invariant as well. Hence all space-time derivatives can enter spin dynamics equations 
only in combinations with the gauge field, i.e. in the form of covariant derivatives: 



As we will see later, this conjecture is correct. 
3 Weak Polarizations 

Spin dynamics equations in the case of weak polarization have been formulated by 
Silin and Leggett [)7| in the framework of the Landau theory [g] . Leggett showed that 
the first two harmonics of the distribution function, viz., the spin and the spin current 
densities, decouple from the rest for slow enough spatial variations. The purpose of 
this section is to generalize the Leggett equations to the case of an electric field present. 

The Landau theory of Fermi liquids operates with phenomenologically defined quan- 
tities. So it is not convenient to literally use the gauge invariance approach discussed 
above. We will set this approach aside until Sec. V and formulate the theory in the 
usual terms of the quasiparticle distribution matrix hk and the quasiparticle energy 
matrix ik- 

In the absence of fields the equilibrium values of hk and Ck are 



where m* is the quasiparticle mass. 

According to Landau's ideas the departure from equilibrium distribution 6hk changes 
the quasiparticle energy along with external fields: 



D.X = d.X + gA.xX. 



(6) 




6ek = de""^' + Y.^p'[fkk'Shk'] 



(8) 



k' 



where (5e|^* is the energy of a quasiparticle in external fields, fkk' = ^F^i^kk'^^' + 
F^i^,&&') is the quasiparticle interaction and Np = m^/cp/TT^/i^ is the density of states 
on the Fermi surface. 
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neia: 

= -fiH = (9) 

As the microscopic Hamiltonian (g) indicates, an electric field alters the equilibrium 
(kinetic) part of the energy matrix by shifting all the momenta of quasiparticles: 

^ + ^e^, where ^e^ = —vpifiAi. (10) 

As a consequence, the equilibrium distribution = fo{ef) will also change: 
fi^ + {dfo/d^)6e^, which, in its turn, will alter the departure from equilibrium Sfik 
defined as Sfik = fik — n\: 

5fik 5fik - -^^^k- (11) 
The amendment in 6hk through relation (g) will change 6ek'- 

Sh ^ ^^fc + E Sp'ihk^ Set,] = 5h + ^5el. (12) 

Now let introduce macroscopic quantities. Spin and spin current densities are de- 
fined according to the formulas: 

= ^ESp[o-nfc], = ^ESp[6-nA:Vfeifc]. (13) 

Both of these expressions in equilibrium in the absence of external field are equal to 
zero. Expanding them to within the first order in departures from equilibrium yields 

s = hj:sp[&6hki (14) 

k 

J I = ^T.^P[^ {^nkVk,el + hlVk,Sek)] = ^Y.VFi^p[o- (^^^ (15) 

where the Fermi velocity vpi = hki/m* and we took the second term in ( p^) by parts. 
This term describes backflow (cf. [^). The expression in parentheses in ([T^) can be 
easily seen to represent the departure from the local equilibrium n|?*^ = /o(e^ + Sik) ~ 

When an electric field is switched on both spin and spin current densities in equi- 
librium will still be zero. This follows from the fact that the amendments to their 
magnitudes brought about by 6h^ and ^e^ are described by eqs. (|T^) and (|T5|) respec- 



tively. But these expressions vanish when Sh^ and ^e^ are substituted for 6hk and 
dsk- Thus an external electric field produces no spin current in equilibrium. We will 
exploit this fact later when writing out the collision integral in the kinetic equation in 
the r-approximation. 

Substituting the expression (g) along with the renewed Sfik ([Tl]) and Sck ([T2|) 
into ([T5|) yields 

J^ = -^{VF. Sp[&5eT'-'^-])f^ + 2 (l + y) ^ Sp[^^n,], (16) 



ratea into an eriective externai-neias energy 



.ext.eff. _ ( . . , /i , ^'t 



^er"" = - Y ( + ^F. ( 1 + ^ 1 gA, 1 (17) 



On plugging this into (|TB|), we obtain the following expression for the spin current: 

Jz = J- + Y)E^Fz^/fc, (18) 

where crk = \ Sp[6-(5nfc]. As compared to the case of no electric field [|7|, the expression 
for the spin current contains an additional term 



^pU XnW^ 



where p = /cp/37r^ is the hquid's density, = ^f(1 + + ^f/3) and Xn = 

g^Ti^ Ny / + Fq) is the normal Fermi- liquid susceptibility. This term appeared pre- 
viously in the spin supercurrent in superfluid Fermi liquid [0] . 

We now proceed to the derivation of the dynamics equations themselves. 

Let expand hk and in the basis of the identity and the Pauli matrices: 

fik = /fci + cTfcO", Ck = efci + efcO". (20) 

In the weak-polarization case, where spin polarization is a very small fraction of the 
total number of spins in the liquid, fk and may be put equal to their equilibrium 
values (|7D (see 0). Then from eqs. (g), (P), ([ig) and (P) we obtain 

efc = --U^L + ^Fz 1 + + — E^fcW- (21) 



2V V 3/" 7 Nf 



k' 



Note that here the electrical-field-dependent corrections to and Sck also appear 
in effect as if it is the external energy ^e|^* which contained them. What we mean is 
that the above expression could be obtained from the expression (§) with the effective 
external energy ([T7|) without allowing for the corrections (|TDD and (|T2|). It is tempting 



therefore to attribute the term in ( fTTp that depends on electric field to an external 
spin-orbital energy. Nevertheless, since it contains Fermi-liquid constants this term is 
not, strictly speaking, an energy in an external field. Quite the reverse, it appeared 
from corrections to the kinetic energy. Thus we conclude that the parallel drawn is 
formal. 

The Silin-Leggett [j7(] kinetic equation for the spin part CTk of fik is (for brevity we 
omit A;-indices) 

dtCT - '^ex CT + VYiVi [cT - = {dtCT)co\\. (22) 

Note the wrong sign before the second term in |7|] . 
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m^k{r,t) = ^jd^cTk{r,t), (23) 



2 

we can further reduce the kinetic equation to 
4 

liNx 

Ff\ ^ TlNy ^ ( ( Fi 



(r \ TliS ( ( r \ \ 

1 + ^1 gA, X J-VF,V^ f a;L + vpj f 1 + ^ j gAA = (^^mjgi) 

Following Leggett, we seek after a solution in the form of the sum of the zeroth and 
first harmonics of m^: 

?im^ = <S' + -^(^l + ^) '^,(J,-jO). (25) 

The higher harmonics can be proved to decouple from the first two if only the char- 
acteristic spatial scale A is greater than the mean free path or the molecular field 
length |7|: 

A > min < -uft, — ^ \ . (26) 
I ^cjl J 

Here k = (Ff / 3 — Fq^) /(I + F^) . This constraint coincides with that for the case of zero 
electric field because the "electric length" Ie = {gE/c)~^ appearing as a new scale in 
the problem is much greater than all the above mentioned scales (for an electric field 
as strong as £" = 3 • lO'^V/cm the "electric length" shortens only to Ie = lO^cm). 

When the mean free path is shorter than the molecular field length (which obviously 
is equivalent to the inequality kcoi^t < 1) we deal with the hydrodynamic regime. The 
opposite case corresponds to the collisionless regime. 

After some algebra we obtain a set of two equations: 

DtS + D,J, = 0, (27) 

2 2 T 

AJ^-5tJ? + ^A(5-5^^) + /s:^5x J, = (28) 

Here S^^ = Xh^l/s'^ is the equilibrium magnetization, t\ = t j (l+Ff/3). As compared 
to the Leggett equations, all the spatial derivatives in (P7|), ( |2S| ) are replaced with the 



covariant ones. The term c^iJ^* = (xn/5'^)('W^/3)(Via;L — VcjlO is non-zero only for an 
electric field varying in time. We will assume this is not the case in what follows. 

First, from now on we make use of the units wherein Xn = S'^- Next, we transform 
to the frame, rotating with the local Larmor frequency u3\Xr) that we assume to be 
parallel to z. To that end we substitute S{i) with lZ{t)S{t) and similarly for J^, 
where IZ = lZ{—uJi^t) is a 3D rotation matrix. Then we get from 



dtS + AJz = 0, (29) 
dtJ^ + —D,{S-u:i^) + i^SxJ, = \ (30) 

3 Tl 



6 



We here used the tact that Caij f-Cif3 f-Cj-y = K-ak^kp-y, which tohows trom the general 
equaUty holding for an arbitrary matrix: eaijT^akT^ifiT^i-i = Cfc/j^ det 

Eq. ( pop is an inhomogeneous linear differential equation on Ji. Its solution can 
be represented as the sum of the general solution of the corresponding homogeneous 
equation and a particular solution of the inhomogeneous one. The general solution of 
the homogeneous equation contains the factor exp(— t/ri) and dies out in a time ri. 
Thus for the most of the time of an experiment we may take into account only the 
particular-solution part of the spin current. 

In the absence of electric fields Di = and both the inhomogeneous term and 
the coefficients in Eq. are independent of time. In this case Eqs. ( |3D|) allow 

exactly stationary solutions, wherein dtS and dtJi are exactly zero in the Larmor 
frame. 

In the presence of an electric field, the inhomogeneous term in ( |3D| ) contains, besides 



the part constant in time, a fast contribution from the rotation matrix that oscillates 
with the Larmor frequency cjl- Hence the solution will be the sum of constant and 
oscillating terms. The oscillating term in the Larmor frame will appear from the point 
of view of the laboratory frame as a doubled frequency motion and must induce a 
doubled frequency harmonic in the NMR signal. 

However the amplitude of this harmonic as well as the correction to the stationary 
motion due to electric field will be of the order of smallness of gEL/c^ where L is the 
spatial scale of the problem. Therefore we can expand S and Ji into degrees of electric 
field to within the first order: S{E) = S + 6S, and similarly for Ji, where S and Ji 
now designate the corresponding quantities in the absence of electric field that meet 
the conventional Leggett equations. For the departures 5S and 6Ji proportional to 
the electric field we obtain from (|2^), ( |3D|) 

dtSS + V,(5 J, + (/A, X J, = 0, (31) 

7/;^ ~ ST- 
dtSJ, + —(ViSS + gA,x(S-u;i^)) + K(SSxJ, + SxSJ,) = -. (32) 

3 Tl 

Thus 5S and 6Ji (deemed as components of one unknown quantity) satisfy a system 
of inhomogeneous linear differential equations with coefficients depending on coordi- 
nates but independent of time. The electric field enters the equations only through 
inhomogeneous terms and hence the corresponding homogeneous system describes the 
time evolution of small perturbations to a known unperturbed zero-electric-field solu- 
tion. The general solution of the homogeneous system is the sum of particular solutions 
in which 6S and 6Ji depend on time through the factors e~*^*. The perturbation fre- 
quencies cu are to be determined by solution of the homogeneous system with the 
corresponding boundary conditions. For the unperturbed distribution to be stable the 
imaginary parts of all the possible frequencies cu need to be negative. Then the arisen 
perturbations will die out exponentially. 

To find the quasistationary part of the solution (i.e. the one that is constant in time 
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Figure 1: Geometry of idealized spin- diffusion experiment. The circles are turns of wire of the rf coil 
used to invert the magnetization in the left reservoir. 

in the Larmor frame), we should average the inhomogeneous terms over the period of 
precession I/cjl, which yields (Aj) = A^z. 

We remind that = —c^^CzikEk- Thus the 2;-component of the electric field 
have no effect on the quasistationary dynamics. If we choose E \\ x for definiteness, 
then Al = c^^Ey^ and the electric field falls out from the equations for the Xi- and Zi- 
orbital components of the spin current, i.e. the electric field affects the quasistationary 
motion only if the ^^-components of the spin current is non-zero. 



4 One-Dimensional Flow 



As an example we consider spin flow through a thin channel directed along the ^j-axis. 
But first we examine the spin distribution in such a geometry in the absence of an 
electric field. The situation is reminiscent of the flow along the z^-directed channel 
studied experimentally [|TD| and theoretically [|n| . We may extend the results obtained 
in Ref. [1TT| to our situation since in the derivation of distributions and in the stability 
analysis in Ref. []rT| only the fact that the flow is one-dimensional is used and the name 
of the coordinate along which the distribution changes makes no difference. 



4.1 In the Absence of Electric Field 

In the idealized geometry of Ref. [|rT| two reservoirs of a Fermi-liquid polarized by an 
external magnetic field directed along z are connected via a long thin tube of length 2L 
and cross-sectional area A <^ L^. After a 7r-pulse of a rf magnetic field, inverting the 
spins, is applied to one of the volumes, diffusion tends to diminish the large longitudinal 
polarization gradient along the channel and eventually equalize polarization in the two 
chambers (see Fig. 

Since the time L^/D to diffuse through the tube is much longer than the time 
A/D to diffuse away from the tube entrances, we may, following [[TT], assume fixed 
boundary conditions S{y = ±L) = i^oi;, where 5*0 is the magnetization density in 
the reservoirs. For convenience we transform to the dimensionless units: s = y/L, 
r = DqI/L^, m = S / Sq, j = LJy/ DqSq, and introduce the parameter v = —kSqTi. 
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iiquiQS ana tnus v is negative. 

Now we rederive the results of []TT| generalized to the case of the parameter ttiq ^ x 



(see below). 

The Leggett equations in the Larmor frame, which are just eqs. (12^1), ( PD|) with 
Di = Vi, for stationary solutions depending on only one spatial coordinate yield 

aj = 0, {d, + vjx)m + j = 0. (33) 

We get 

j = const, m(s) = lZ{—svj)ino — sj, (34) 

where mg = m(s = 0) is the magnetization in the center of the channel. Its absolute 
value is \/l — j'^ and its direction remains an arbitrary parameter. This uncertainty 
in the direction of rho corresponds to the indefiniteness of the initial phase of the 
precession. Actually the result in [tTT[] is written out (with some slips) in the coordinate 
representation corresponding to the choice of mg || x. This is why we dwelt on the 
derivation. 

Then the boundary conditions rnl^^^i ~ require lZ{2vj)mo = —mo. Hence 
either mo = and we get 

j = —z, m{s) = sz. (35) 
Or j _L mo 7^ and vj = 7r/2 + Trn, n G Z. The stability analysis and numerical 



simulations [JTTj] show that only the solutions with n = and n = —1 are stable. For 



liquids with v negative n = — 1 is relevant, while n = corresponds to v positive. Using 
the overt form of the rotation matrix ^{0) = {6af3 — 0^0 13) cos 6 + 6a0p — e^^^^^sin^, 
we obtain for n = — 1: 

j = — — , J = —jz+niQXz, m(s) = — moXjsin(7rs/2)+mo cos(7rs/2)— sjj. (36) 

In such a solution the polar angle tt — 'd between j and z is fixed by the condition: 
j = —{jz) = COS'S. (Here we see that of necessity j < 1 and \v\ > 7r/2.) The arbitrary 
azimuthal angle of j is parameterized by a constant unit vector mo = z x j/ttiq. 

The three unit vectors mo x ''^o, 3 form a right-handed basis (see Fig. 0). 

Subject to the value of \v\ either of the above stable stationary solutions takes 
place. For \v\ < 7r/2 the solution is the longitudinal diffusion (|35|). For l^;] > 7r/2 this 
configuration becomes unstable (Castaing instability) and ( |3^ ) is realized. 

Properly speaking, arbitrary boundary conditions can be considered in the same 
manner but only the ones involved, i.e. of the diametrically opposed magnetizations 
in the two reservoirs, appear to have a solution of more or less simple form. 

In fact, the diffusion under consideration is only ^wasz stationary. That is, it is 
independent of time on the short scale LP' / Dq. But on long times owing to the finite 
size of the reservoirs, slowly as it will, the current will change the magnetization 
densities in the reservoirs, thus changing the boundary conditions. Thanks to this 
slowness we may deem the situation at any given moment as stationary on times ~ r 
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Figure 2: Disposition of the right-handed basis in the non- longitudinal one-dimensional stable stationary 
solution of Leggett's equations with the boundary conditions as on Fig. 1. This solution is reahzed 
when \v\ > tt/2. The unit vector mo is directed from the plane of the picture (the j — z plane) towards 
the reader. The angle i) equals i) = arccos(7r/2|f |). It varies from -(9(0) G (0, vr/2) at the start to zero at 
the time of switching to the longitudinal solution t = t* — T/2. 



and simply substitute Sq for z in the boundary conditions, so that now they read 
m{s = ±1) = ±So. 

For the long-time dependence of the magnetization in the right reservoir we then 
get 

Tfi = 5o,-, (37) 

where T = VL/DqA is the long-time scale. Thus the quasistationary approximation 
is valid as long as T ;> L^/Dq, or, equivalently, V ^ LA. For experimental conditions 
of |D| T - 100 sec. 

For l-ul < 7r/2 the current j = —Sq and the magnetization decays exponentially: 



SoW = So(0)exp(-f/T). 



(38) 



For \v\ > 7r/2 the magnetization density diminishes algebraically until it reaches 
the critical value 7r/2. It is convenient to introduce a constant ( = 7t/2kti so that 
j = C/Sq. Upon substituting j = j{—jSo + \/l — prriQ x Sq) eq. ( |37| ) yields 



dt Sq 



T 



dSc 

dt 



1 

^0^ 



1 - ^rno X Sq. 

Dq 



(39) 



Then while v{t) exceeds the critical value, i.e. for times t < t* — T/2, where 
r = TSl{0)/2e. we get 



So{t) = 



t* - t 



\ T 



SQ{t) = n 



( 



niQ 



1 
Sq\ 



1 



C_dt 

<2 




ni7fiQpit))z, (40) 
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Figure 3: Disposition of the right-handed basis: a) for the non- longitudinal solution at arbitrary non- 
zero times and b) for the longitudinal one following the non-longitudinal. The angle (3 becomes non-zero 
during diffusion, at the start Sq = z. 



where 




t 



1 



J 



(41) 



The absolute value of magnetization decays as a square root of time and its direction 
rotates around rho, i.e. in the j — z plane (see Fig. (§a)). The angle (3 increases 
monotonically with time. As numerical evaluation shows, it reaches tt at the critical 
point for the first time when the duration of the non-longitudinal part of diffusion 
constitutes t* « lOT. 

Note that the time evolution of the magnetizations in the two reservoirs is symmetric 
as a result of the conservation of current. And at any moment the directions of the 
magnetizations in the vessels remain opposite to ensure the diametrically opposed 
conditions involved. 

For t > t* — T/2 the magnetization will proceed by exponential slump (see Fig. (H)) 
with the angle /3f = - T/2) 7^ between z and (see Fig. (gb)). 

4.2 In the Presence of an Electric Field 

Let now an electric field E = Ex be applied to the liquid in the channel at some 
time Ie > In accordance with what was said above, all quantities will acquaint 
increments proportional to the electric field: the magnetization and current in the 
channel m ^ m + 5m^ 3^3+ ^3- And the magnetization in the right reservoir 

^ Sq + SSq. 

The boundary condition for 6m will be Sm{s = ±1) = ±^S'o/*5'o. 

Since the quasistationary inhomogeneous terms {g{Ai)x) are independent of time, 
we seek for a quasistationary solution dr6m = 0, dr6j = 0. The oscillating terms will 
be discussed later. Of course, Sq, m and j all depend on the "long" time t/T. 



11 




Figure 4: The time dependence of the absolute value of the magnetization density in the reservoirs 
So{t) in case the initial value 5*0 (t = 0) exceeds the critical value 7r/2 at which the longitudinal diffusion 
becomes unstable. At t = t* — T/2 the value of 5*0 reaches 7r/2 and from this point diffusion switches 
to the longitudinal regime. 

After introducing the additional parameter e = gEL/c we get 

djj + ezxj = 0, (42) 
{ds + vjx)6m + ez X m + 6j + v6j x m = 0. (43) 

The first of the equations yields 

S3{s) = Sj,{s) + Sjo. (44) 

where 

^Je(s) = X js (45) 

and SJq is a constant to be determined from eq. (^) with the boundary conditions 
imposed. We can make a substitution 6m{s) = lZ{—vsj)6rh{s) into eq. (^), after 
which it can be easily integrated 

s 6 S 
5m{s) = lZ{—vsj) / lZ{vsj){ez x m + Sj + v6j x m)ds -— . (46) 

And then the boundary conditions reduce to 

1 S S 

I 'TZ{vsj){ez X m + 6j + v6j x m)ds = 2Ti{vj)—^. (47) 

As a matter of fact, SJq should be determined from this equality. At the moment 
of switching on of the electric field SSq = and the above equation can be solved 
analytically. SJq then turns out to be zero. Nonetheless, although SJq is, generally 
speaking, non-zero at arbitrary times, this part of 6j is constant in space and does not 
change the symmetric picture of the diffusion. 
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interest to us. inis part nas opposite signs at tne two tube entrances, tnus cnangmg 
the magnetizations in the two reservoirs by the same amount, i.e. not symmetrically. 

Unless j is parallel to the term is directed along mo perpendicular to the j — 2; 
plane. Therefore, as diffusion goes on, the non-symmetric change of magnetization in 
the reservoirs will lead to an additional phase shift between the two vessels. 

To find the change in the time evolution of Sq due to electric field we substitute 
5*0 with Sq + 6S0, where So{t), as previously, belongs to the case of the absence of 
electric field and the part SSq is proportional to the electric field. This part obeys an 
equation, which can be obtained by taking variations of eq. ( |37| ) using 



= So [6jo + Sj,{s = 1)] + 6Soj. (48) 

Next we expand 6S0 onto the components in and perpendicular to the j — z plane: 
^5*0 = 6S]^ + 6Sq^^. For the component 6S]^ perpendicular to the j — z plane we get 
TddSfP/dt = Soldjip + 6je{l)]. The phase shift is described by the term 6je. Hence the 
phase shift increases with time as 



'SQ{t)smP{t) So{t) sin f3{t) 



where we have neglected the presumably small component SSq~^' compared to Sq in 
the denominator. 

For the non-longitudinal stable solution (|3§) we have (5je(l) = ejsm[P(t) — 'd{t)]. 
and the corresponding integral can be found only numerically. For times close to ts 
we have {t > ts) 

t-tE ^sin[(3{tE) - ^{tE] 



1 sm/j[tE) 
When tE <^T^ the last expression simplifies 

5a « -2^—^. (51) 

However, this formula cannot be applied to tE extremely close to zero. Namely, in the 
case that the transverse component of the magnetization in the reservoir SQ{tE) sin f3{tE) 
is of the order of magnitude of 6Sq^, the up-to-the-first-order extension of the transverse 
component of Sq in electric field is not valid. 

For the longitudinal stable solution (|35| ) eq. (^) can be easily integrated and yields 



sin Pf in both the numerator and the denominator, where = P{t* — T/2) here is the 
polar angle between Sq and z handed down from the first (non-longitudinal) part of 
the diffusion. Thus 

6a{t) = -2e[exp{{t - to) /T) - 1] (52) 

for all values of /?/ = [3{t* — T/2) except for [3j = 7rn, n G Z, when the phase shift 
is undefined. In particular, if there were no non-longitudinal regime at all (due to 
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woma De no aaaitionai pnase snirt. 

The time in (1^) is either t* — T when the longitudinal solution follows the non- 



longitudinal regime, or Ie if an electric field is turned on already at the longitudinal 
part of the diffusion. In the former case, eq. ( |5^ ) should be added up with the phase 
shift accumulated during the non-longitudinal part. 

Exp. (^2|) formally tends to infinity with time, but owing to the small coefhcient e 
it still remains very small for all of the time of an experiment. 

The appearance of an additional phase shift in an electric field is quite understand- 
able. The diffusion process may be outlined as a transfer of down spins from the left 
reservoir to the right one. Naturally, the current along of down spins is equivalent 
to the current of up spins in the opposite direction. In accordance with the boundary 
conditions the spins of ^He atoms at the left end of the tube are directed preferably 
down. The "surplus" of down spins drifts to the right and interacts with the electric 
field E = Ex as if there was an increase to the magnetic field: —{EvF/c)fiz. 

On the contrary, the up spins at the right end float to the left, and their spin-orbital 
energy [Ev-p / c) fiz amounts to an effective decrease of the magnetic field. The spins 
at the opposite ends thus precess at slightly different frequencies, unless, of course, 
they are not confined to the z direction like in the purely longitudinal solution. Hence 
in the second solution (for \v\ > tt/2) there develops a phase shift between the two 
entrances of the tube. 

The conclusion about the existence of an additional phase shift was drawn for a 
situation with diametrically opposed magnetization directions at the tube ends as 
the boundary condition. But the inference in itself remains valid for other boundary 
conditions as well, irrespective of the fact that we are ignorant of zero-electric-field 
distribution in that case and of exact expression for the phase shift. This follows, first, 
from the current conservation in the absence of electric field, which secures the self- 
similar boundary conditions in the quasistationary diffusion process. And secondly, 
from eq. (0), which holds for all stationary solutions in one-dimensional geometry and 
leads to non-symmetric change of the magnetization in the reservoirs when an electric 
field is applied. 

As for the oscillating part of the gauge field, it does not contribute to the current 
along the channel. Indeed, the oscillating part of the gauge field is 

gAi = {g E /c)zi{x sin iOi,t — if cosui^t) (53) 

if the electric field is parallel to Xi as before, that is, gAi is proportional to Zi. Taking 
into account that different i-components of 6Ji do not enter each other's evolution 
equation (^2]), it means that in the particular solution of the inhomogeneous system 
only 6S and 6Jz will be non-zero and dJ^ and 6Jy will figure only in the general 
solution of the homogeneous equation. According to what was said above, if only the 
unperturbed solution is stable, the latter dies out exponentially with the rate (Imo;)"^, 
where uj is an eigenvalue of (0),(^)- 
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a line structure across tne cnannei, wnicn vary on a spatial scale oi ana oscillates 
with the Larmor frequency in the Larmor frame. Though in general it is rather com- 
pHcated, the main features of such a structure may be seen on the easier example of 
the longitudinal diffusion, which is considered in the Appendix. 



5 Strongly-Polarized Fermi- Liquid 



Landau theory cannot be literally applied to strongly spin-polarized Fermi liquids [|T2 



However for T = Fomin derived microscopically transverse spin dynamics equa- 
tions fO], the form of which, in fact, coincides with the collisionless limit of the 



Leggett equations. 

In this section we will consider the appHcation of the SU(2) gauge invariance directly 
to a microscopic Hamiltonian of a liquid at T = 0. The Hamiltonian of the liquid 
involved in the absence of external fields consists of the kinetic energy]^ 

Hk = I ^+(r)V2^(r) (54) 



2m 

and an interaction Ti^^i between particles. We neglect the very small spin-orbital 
interaction between particles. 

To include electric and magnetic fields into consideration we must use (g) instead 
of ?-^K- Such a replacement can be interpreted as a result of a gauge transformation 
djj^ = — ig^A^. It is convenient to write the Hamiltonian (g) in the form: 

5n = - j gA^{r)S{r) Sr - j gA,{r)J,{-gA,, r) d'r, (55) 
where the spin and the spin-current operators are defined as follows 

S{r) = ^+(r)|v^(r), (56) 



'iV,ip+{r)) ^il>{r) + (-'iV,'i/'(»-)) + -'i/'+(r)i/'(r)X, 



;57) 



2m 

Following Fomin [|T2| we now transform in each space point to a frame precessing 



at such a frequency 17, that the magnetization in the rotating frame is equilibrium. 
fl is supposed to vary slowly enough in space and time so that we could leave only 
first terms in lj and /c, where cu and k are respectively the characteristic frequency and 
wave number of Q, in the Larmor reference system. 

As a result we get the following expression for the Hamiltonian in the rotating frame 

n = nK + Hint + 6n + J S{r)Ct d^r. (58) 

Now we introduce a and P as the spherical coordinates of — ft: ft = — ^(sin f] cos a, sin P sin > 
and rotate the spin space i/j TZi/j, where 

7^ = exp{-iaay2) exp(-z/3(jV2) exp{-i-fay2) (59) 



^We leave out spin indices of the field operators of '^He atoms for brevity. 
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eter. ine nanay constraint on it wiii be clear later, men tne space-time aerivatives 
d^ip transform as follows: 



(60) 



with the covariant derivative D^il) = d^il) + {7l^d^j7l)ilj = d^il) — i^-uj^ip. This must 
be regarded as the definition of u;^. It can be written out equivalently as 



^ — sin cos 7 d^a + sin 7 d^P ^ 

sin P sin 7 d^a + cos 7 d^P 
^ cosPd^a + 5^7 J 



(61) 



where IZai = lZz{—j)'R^y{— P)lZz{—a) is a 3D rotation matrix. 
Under this transformation 

= Hk + ^int + / S{r)nn - /(a;o + a;L7e.^)S'(r) d'r - 

[uj, + gnA,)J,{-uj, - gUA, r) Sr. (62) 

As can be easily seen, l^z = (— sin/3cos7,sin/3sin7,coS/5), l^Q. = — Qi;. 

In case that the absolute value of 17 is greater than its variations in the reference 
system that precesses at the Larmor frequency (i.e. greater than a + cjl, 7 and 
Vfa/m, V'-P/m, Vf7/m), then the first three terms in ( |52|) are the greatest. Together 
they constitute the unperturbed Hamiltonian of a Fermi liquid in a uniform field —Qz. 
The ground state of this Hamiltonian is {S) = (0, 0, S) and ( J^) = 0. The former is 
extremely important because shows that a and P serve as spherical coordinates of S 
in the laboratory frame. 

The other terms in (^) are the perturbation. Its magnitude in the ground state 
equals 

- (u;o + uJi.nz){S) = -5(7 + cosP{a + c^l))- (63) 

The expression above contains only time derivatives. To include space derivatives 
one should allow for the second order corrections, which can be evaluated as in Maki's 



paper 



AF 



X" 



X 



J 



■[(Va)2 + (V/3)2 + (V7)2 + 2cos/3VaV7] 



^ + X ^aik^c 



^ 2 



Note the wrong sign in Ref. [O]. A new definition is 



Cbi = 71. ^u;,- = 



^ sin /3 cos Of Vi7 — sinaViP^ 

sin/3sinaVi7 + cos aViP 
^ cos/3Vj7 + Via J 



(65) 
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witn respect to 7, we oDtam 

S + x'^V( V7 + COS pVa + S X ^E) = 0. (66) 

In the end we could set a constraint on 7 guided by the reasons of rationality, i.e. so 
that S = 0, which requires V7 + cos /^Va = —S x gE/c. Substituting this condition 
into the equations obtained by varying Ces with respect to a and we obtain a set 
of two equations for a and P which are equivalent to 

DtS = ^D,{D,SxS). (67) 

It is interesting to trace the formal transition from Eqs. ( |27]) , (|2B1) to Eq. ( |57|) 
though the modified Leggett equations, definitely, cannot be applied to a strongly- 
polarized liquid. 

First of all, for a strongly-polarized liquid S is much greater than cjl, hence, S^^ 
could be neglected in (|28|). Moreover, the collision integral Ji/ri vanishes as compared 
to the molecular-field term kS x by virtue of the parameter kSti ^ 1. Thus we 
can get to the limit kSti 1 in an expression for the quasistationary current, which 
can be obtained by resolving Eq. ( |2S| ) with respect to Ji provided that DtJi = 0: 

J^ = - i^^^snf "^^'^ ^ ^^^^ ^ {^TifS{SV,S)]. (68) 

As was discussed concerning the weak-polarization case, the transition to the qua- 
sistationary current is self-consistent only for the quasistationary part of the solution 
because only then both sides of (|681) are constant in time in the Larmor frame. The 
senior (the third) term in ( |65| ) identically vanishes for homogeneous spatial distribu- 
tions of the absolute value of magnetization. We assume this condition to fulfil, for 
the collisionless regime of the weak-polarization case it can be proved. 

Substituting then the expression for the quasistationary current into ( P7|) , we ob- 
tain ( |57|) with xj = w'^/Sk, or, in ordinary units, xj = {Xn/ g^){w^ /"^i^)- It is this value 
of Xj that the Fermi-liquid theory gives in the limit of weak polarizations [|T2|. Still 
Xj cannot be obtained in the general case. 

From the fact that the modified Leggett equations (|27|), ( |2S| ) yield (|571 ) in the deep 
collisionless limit in the quasistationary approximation it follows immediately that 
all conclusions about the quasistationary part of spin flow through a channel and, in 
particular, of the additional phase shift remain valid in the strong-polarization case. 



6 Conclusion 

In this paper we considered how an external electric field can influence spin dynamics 
of an electrically neutral Fermi liquid through the spin-orbital interaction with nuclear 
magnetic moments. In the framework of Landau's theory of Fermi liquids attribut- 
ing an additional energy due to spin-orbital interaction to each quasiparticle leads 
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requirements imposea Dy intrinsic oU(^zj gauge mvariance or tne interaction or a mag- 
netic moment with electromagnetic field. The corrections caused by electric field in the 
Leggett's equations are responsible for an additional phase shift in a one-dimensional 
spin fiow experiment. This phase shift, which is proportional to the electric field, co- 
incides in the order of magnitude with that in superfiuid ^He-5 but grows with time. 
Its magnitude is observable experimentally. 

For strongly spin-polarized Fermi liquids, electric field is incorporated into Fomin's 
equations for transverse spin dynamics of a strongly spin-polarized Fermi liquid at zero 
temperature. The formal correspondence between these equations and the coUisionless 
limit of Leggett's equations spreads to the case of an electric field present as well. 
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Appendix 

In the case of longitudinal diffusion the inhomogeneous term in ( pT]) is identically zero, 
but not in (|32|). Thus the solution of the coupled system ([1?T]) , ( |32| ) should as usual be 



sought in the form of the sum of the two terms depending on time through the factors 
sincjL^ and coscjl^- 

According to what was said in the body of the paper, we seek for a solution of the 
system 

dt6S + d,J, = 0, (69) 

dt6J, + — (d,S^—( X sin ULt-y cos ui,t)x(S-u>L)) ^k6Sx J, = -^0) 
3 c Ti 

in the form 6Jz = {DqSq/ L) (a sin uji,t + b cos uji,t). From the first equation we get 
6S = {DqSq/ L){—{hl / uj\)j sincjL^ + [a' /uj\,) coscjl^), where a' = dza. 

Next we should utilize the boundary condition of the absence of spin transfer 
through the walls of the channel: 5J z = at the walls. Hence the solution can 
be expanded into harmonic modes, in which 5Jz depend on z through the factors of 
the form sinkz. Here k = irZ/d and d is the thickness of the channel, which in general 
is a function of x. 
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a = 



e(s - ujl/Sq) 
C2 + 1 



{Cx + y) 



b = 



e{s - ujl/Sq) 
C2 + 1 



(x - Cy) 



(71) 



where C = cjlTi — Dq^^/'^l — vs. We here allowed for the fact that due to diffusion 
after some finite time after the beginning of the experiment is not equal to ui,. 
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